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ABSTRACT 

This paper presents recent work conducted on the integration of 
mass-interaction physical models in the Faust programming lan¬ 
guage. After a brief introduction to mass-interaction networks, 
FAUST, and previous works on this topic, we present a simple mod¬ 
eling framework, a FAUST code generator and its associated library, 
allowing to implement ID mass-interaction models. In addition to 
the open-source tool itself, this research offers a perspective on for¬ 
malizing arbitrarily large networks of bidirectional feedback cou¬ 
plings and state-space models in FAUST, through routing patterns. 
We finish with a set of examples, and discuss future perspectives and 
challenges. 


1. INTRODUCTION 

For several decades, physical modeling has been used to synthesize 
audio by means of simulating the behaviour of vibrating objects. A 
panoply of methods have been proposed over the years, from lumped 
discrete models m. to Waveguides [2|, to large scale Finite Differ¬ 
ence schemes m, that have gained in popularity with the increase 
of computing power. Creating a model of a mechanical instrumental 
system can be simpler than explicitly formulating the signal that it 
produces (as sound properties emerge from the physical conditions 
of the matter) and offers direct means for control and interaction, ei¬ 
ther by simulating musical gestures or by coupling the user and the 
virtual object, for instance using haptic technologies |4l . 

Faust 0 is a functional programming language for real-time 
Digital Signal Processing (DSP) with a strong focus on the design of 
synthesizers, musical instruments, audio effects, etc. The Faust 
compiler can be used to “translate” a Faust program to various 
non-domain-specific-languages such as C++, C, JAVA, JavaScript, 
LLVM bit code, WebAssembly, etc. Thanks to a wrapping system, 
code generated by Faust can be easily compiled into a variety of ob¬ 
jects ranging from audio plug-ins to standalone applications, smart¬ 
phone apps, web apps, etc[j This mechanism also makes it possible 
to add MIDI, OSC, polyphony, etc. support to any FAUST-generated 
program. 

1.1. Mass-Interaction Physical Models 

Pioneered in artistic applications by the CORDIS-ANIMA system 
(Tj at ACROE, mass-interaction physical modeling allows to formu¬ 
late physical systems in the form of lumped networks, composed 
of two main components: masses, representing material points in a 
given space (ID, 2D, 3D) with a given inertial behaviour, and inter¬ 
actions, each representing a specific type of physical coupling (i.e., 

'The FAUST website contains an exhaustive list of all the FAUST targets: 

https://faust.grame.fr 


visco-elastic, collision, non-linear, etc.) between two mass elements. 
Mass-interaction systems are now used in a variety of contexts (mu¬ 
sical & other), partly for the fact that arbitrarily complex virtual ob¬ 
jects can be described simply as a construction of elementary physi¬ 
cal components. A basic model is shown in Figure 1. 
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Figure 1: Topological representation of a mass-interaction model. 
Here, a fixed point (represented on the left) is connected to a trian¬ 
gle composed of masses and dampened springs. An input module 
interacts with the top mass through a non-linear pluck interaction. 


Unlike FDTD methods (3), creating physical models with this 
formalism avoids the need to explicitly define a mathematical model 
(partial difference equations systems, boundary conditions, etc.) for 
a given physical structure beforehand. Therefore, it lends itself par¬ 
ticularly well to iterative and exploratory design of "physically plau¬ 
sible" virtual objects, grounded in the laws of Newtonian physics 
but not necessarily limited to the mechanical constraints of the real 
world. 

Mass-interaction physical models can contain anything from a 
couple of physical elements to tens or hundreds of thousands of 
them. Assembling and configuring the models element by element 
can be very time consuming. To this end, user-friendly modeling 
environments have been proposed, namely GENESIS 0 (and more 
recently Synth-A-Modeler Q) for ID audio applications. The for¬ 
mer offers high level tools for generating topological structures, and 
analyzing/tuning physical constructions through modal analysis 0. 

1.2. Current State of Physical Modeling in Faust 

Various projects have been using FAUST to implement physical mod¬ 
els of musical instruments. 

The Faust-STK 0 is a complete re-implementation of the 
waveguide and modal models of the Synthesis ToolKit (STK) (sa¬ 
lt also contains various models from the Soundius Project. [^] 

2 Unfortunately, there is no documentation/publication on this project yet. 
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Julius O. Smith implemented a series of waveguide meshes that 
landed in the FAUST librarie^Jbut that were never documented/pub¬ 
lished. 

More recently, the FAUST physical modeling toolkit m was 
introduced. It is based on a library allowing for the implementation 
of bi-directional block diagrams in Faust and containing a wide 
range of musical instrument parts that can be assembled in a mod¬ 
ular way. It also comes with mesh2faust tm a tool to generate 
modal physical models compatible with the Faust physical model¬ 
ing library using Finite Element Analysis (FEA). 

The work presented in this paper was partly inspired by Ed 
Berdahl’s Synth-A-Modeler (7] (which itself direcly draws upon 
CORDIS-ANIMA Q] and GENESIS |6j). This environment allows 
for the implementation of hybrid models combining mass-interaction 
systems with waveguide models using a graphical user interface 
(GUI). Synth-A-Modeler is based on a series of FAUST libraries and 
generates custom FAUST programs corresponding to the models im¬ 
plemented in the GUI. While it successfully combines various types 
of modeling techniques at a high level and facilitates their control 
using custom haptic interfaces such as the FireFader ca. it has, 
to our knowledge, never been used to implement large scale mass- 
interaction models. 

Our proposed approach does not aim to supplant Berdahl’s; 
rather, from a similar starting point it questions how the FAUST lan¬ 
guage’s versatility can be used to formalize arbitrarily large mass- 
interaction models - and more generally speaking complex feedback 
networks - in a direct, concise and clear manner. 


Equation 0 can be normalized to unity, and rearranged in order 
to express the mass’ position update scheme (discrete-time positions 
and forces are noted X and F ) : 

X (n+1) =2X (n) -X (n _ 1) + F ^ (3) 

With M, the discrete time inertial parameter defined as : 


Hence, the basic discrete-time mass module produces new po¬ 
sition data based on its current position, previous position, the 
"discrete-time" mass parameter M, and the sum of forces applied 
to the mass from the previous interaction computation step. 

The initial position X(o), delayed initial position X(_i) (which 
infers initial velocity) and initial force T) 0 ) must be supplied at the 
start of the computation. 

2.1.2. Discrete-time implementation of a dampened spring 

The elastic force applied by a linear spring with a stiffness k and a 
resting length of (o = 0 connecting a mass m2 at the position x 2 to 
a mass ml at the position xi is given by Hookes law: 

f s i ^2 = ~k.(x 2 - Xi) (5) 

The exact equivalent of this equation in discrete time is : 


2. MASS INTERACTION PARADIGM IN FAUST 


T’sl—^2(n) — F.(X 2 (n) -^l(n)) 


( 6 ) 


Before getting into implementation specifities, this section presents 
the basics of mass-interaction networks, in the case of ID systems, 
in which all masses vibrate along a single z axis. These models are 
sometimes referred to as "zero-D", since they are purely topological 
and contain no direct geometrical information. First, discrete-time 
mass and interaction physical algorithms are presented and assem¬ 
bled into an explicit computational scheme. 

Then, relying on a matrix-based representation of the topological 
network, we present a generic FAUST architecture that implements 
this computational scheme. 

2.1. Discrete-Time Physical Algorithms 

Below, we present finite difference implementations of two of the 
most basic elements in a mass-interaction network: punctual masses 
and springs. 


Where the discrete-time stiffness parameter K = k. The fric¬ 
tion force applied by a linear damper with a damping parameter z 
connecting the same two masses is : 


/dl-i-2 


d(x 2 —xf) 

dt 


(7) 


Using the Backward Euler difference scheme, the frictional force 
can be formulated in discrete-time as : 


/dl—>2 (f) = -Z. 


{x 2 (t)—x\(t)) - (x 2 (t-AT ) 

— 


— xi(t— AT)) 

( 8 ) 


Which after normalization becomes : 

Fd 1—>2(n) = — ^.((X 2 ( n ) — A 2 ( n _l)) — (A'i( n ) — A' 1 ( n _ 1 ))) (9) 


2.1.1. Discrete-time implementation of a punctual mass 

The motion equation for a continuous time mass is given by New¬ 
ton's second law: 

f = m.a = (1) 

J dt 2 

Where / is the force applied to the mass, m is its inertia a its 
acceleration, and x its position. Applying the second-order central 
difference scheme, with the sampling interval noted AT, a discrete 
equation of the mass can be formulated as follows: 


/(f) = m. 


x(t+AT) — 2x(t) + x(t— AT) 
AT2 


( 2 ) 


With the discrete time inertial parameter Z defined as : 


X = 


AT 


( 10 ) 


The global equation of the force applied by the dampened spring 
is composed of F s and Fd : 


F( n ) — K.(X 2 („) Xi( n )) 

- Z.{(X 2(n) -X 2(n _ 1) )-(X 1(n) -X Hn _ 1) )) 

It is applied symmetrically to each mass (Newton's third law): 

F 2 — >l(n) = F(n) 

-Pl->2(n) = +T( n) 


"https://github.com/grame-cncm/faustlibraries 


( 12 ) 
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2.1.3. Discrete mass - dampened spring - fixed point oscillator 

A linear harmonic oscillator is obtained by combining equations 0 
and dill, in the case where X\ is a fixed point set to X\^ = 0, 
n G ’£. This results in : 

x <"«> - ( 2 -t)-*<"> + + T (13> 

Since the basic oscillator is a very common element in modeling, 
the integrated form given in (JT3J can be implemented in the form a 
specific mass-type module (although it is identical to assembling a 
mass, dampened spring and a fixed point). 

2.1.4. Generalization 

Any element in a mass-interaction model follows the basic template 
of the elements described above. More complex interactions stem 
from conditional statements (e.g. springs only active during inter¬ 
penetration of two material points, as in visco-elastic collisions) or 
dynamic stiffness or damping parameters that depend on the posi¬ 
tion and/or velocity of the connected material points (e.g. through 
non-linear lookup tables, such as in plucking or bowing interactions 

El). 

It is important to note that the M and Z parameters are depen¬ 
dent on the sampling interval. Hence, the oscillatory behaviour of 
physical models will be dependent on the sampling rate of the simu¬ 
lation. 



¥ 


Figure 2: Computation cycles of the model presented in Figure [ 7 ] 
At each time step, the mass-type algorithms are first computed using 
the forces calculated in the previous step, then the interaction-type 
algorithms are computed using the new positions. 


2 . 2 . Computation Scheme 

Computing a mass-interaction model consists in calculating the 
mass-type and interaction-type algorithms in a closed loop. The ex¬ 
plicit time step increment is carried by the masses, as shown in the 
discrete-time equation 0- The interactions in themselves are delay¬ 
less operations, but can be computed since their output is fed back 
into the masses for the next calculation step (cf. Figure [2). In other 
words, calculating a step of real-time audio requires to run all the 
masses’ algorithms once, then all the interactions’ algorithms. 

2.3. Representing the Topological Network 

The topological connections of a mass-interaction model can be for¬ 
malized as a routing matrix of dimensions J x 2 K, where J is the 
number of material elements (or M points) in the network, and K is 
the number of interactions (each interaction module has two connec¬ 
tions -or L points in the usual terminology ID) : 

*0J1 *0J2 ■ • ■ ikjl 

mi / 1 or 0 1 or 0 . 

m 2 1 or 0 . 

mj \lorO . 1 or 0 / 


ik_l 2 

1 or0\ 

(14) 


As an example, 0 presents the routing matrix for the topolog¬ 
ical structure shown in Figure [T] The material elements (fixed point, 
three masses and a position input module) are represented vertically 
and the L points of the four springs and the non-linear interaction are 
represented horizontally. 

The closed-loop physical calculation scheme performed by 
FAUST is shown in Fig. [3] On the left, a LinkToMass connec¬ 
tion function routes the force feedback signals produced by the in¬ 
teractions based on the routing matrix (thus calculating the sum of 
forces for each mass). The new positions of the material elements 
modules are then calculated. These positions are then fed into a 
MassToLink connection function, that routes the signals to all of 
the concerned interactions. Finally, the pairs of force signals pro¬ 
duced by the interactions are fed back for the next calculation step. 

Position and force inputs are directly incorporated into the 
LinkToMass function, so that they are applied to the correct in¬ 
put module. Similarly, modules whose positions are observed as 
audio outputs are simply added as extra signals at the end of the 
MassToLink function. 


2.4. Faust Implementation of Mass and Interaction Elements 


Each column in the matrix must have a single connection set to 
1 and all others to 0, as an L point only connects to a single M point 
(partially connected interactions are not allowed). On the other hand, 
a material point could be connected to any number of interactions in 
a given model (many connections set to 1 for a single line). 


The mi . lib library contains the Faust implementation of most 
elementary mass-type elements (i.e., masses, fixed points, oscilla¬ 
tors, etc.) and link-type elements (i.e., springs, collisions, non-linear 
plucking / bowing, etc.). Since the implementations are similar, we 
will explicit only the two simplest and most common elements be¬ 
low: the mass and the spring. 
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Figure 3: Faust -generated diagram corresponding to the model presented in Fig. [ 7 ] 


2.4.1. Mass 

The discrete-time algorithm of a basic mass module described in 
0 can be easily expressed with letrec environment expression 
in Faust: 

mass(m,xO,xl) = equation 
with { 

A = 2; 

B = -1; 

C = 1/in¬ 
equation = x 
letrec{ 

'x = A* (x : initState(xO)) + 

B* (x' : iniState((xO,xl))) + 

* (C) ; 

}; 

}; 

Listing 1: The discrete-time mass algorithm in FAUST. 

The module takes an input signal (the sum of all forces fed back 
through the interaction feedback loop and routing function) and pro¬ 
duces a position output. The initial position and delayed position of 
the module are dealt with using the initState function, which 
initializes the first step with the correct values. 

2.4.2. Dampened spring 

Similarly, a visco-elastic spring expressed in Faust is shown in List¬ 
ing [2] Interaction modules such as the spring take two input signals 
(the positions of the masses connected together by the link) and pro¬ 
duce two identical and opposite force signals. 

Attention must be paid to the correct initialization of velocity 
based interactions, especially when the initial position or speed of 
the masses is non-zero. To this end, the delayed initial positions 


of the two connected mass elements are supplied as arguments to the 
interaction function, which initializes them with the initState ( ) 
function. 

spring(k, z, xlrO,x2r0,xl,x2) = 
k*(xl-x2) + 

z* ( 

(xl - (xl' : initState (xlrO))) - 
(x2 - (x2' : initState(x2r0))) 

) 

<: *(- 1 ) 

Listing 2: The discrete-time dampened spring algorithm in FAUST. 


3. CREATING MODELS WITH MIMS 

Mass Interaction Model ScRiPTEsQis a simple graphical or 
command-line tool written in Python to generate structured Faust 
code from a textual description of a physical model. 

Models are described in a format similar to the PNSL language 
GUI: each physical element has a specific label, specific physical 
parameters and/or initial conditions, etc. Parameters can be added 
to this description and shared by any number of physical modules, 
allowing global variation of the physical attributes (i.e., stiffness, 
damping, mass, etc.) of a subset of modules in real-time. 

MIMS’ physics2faust tool compiles the model by : 

• parsing all of the physical modules and noting any specific 
elements (i.e., position or force inputs, audio outputs, etc.) 

• creating the routing matrix and translating it into the two dual 
FAUST routing functions. 

4 https://github.com/mi-creative/MIMS 
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• ordering the resulting data into the output .dsp file. "Place¬ 
holder" functions are created for position / force inputs, al¬ 
lowing the user to describe his input functions directly in the 
Faust code. 

# Define global parameter attributes 
@m_K param 0.1 

@m_Z param 0.001 

SnlK param 0.05 
SnlScale param 0.01 

# Create material points 
@m_s0 ground 0. 

@m_m0 mass 1. 0. 0. 

@m_ml mass 1. 0. 0. 

@m_m2 mass 1. 0. 0. 

# Create and connect interaction modules 
@m_r0 spring @m_s0 @m_m0 0.05 0.01 
@m_rl spring @m_m0 @m_ml m_K m_Z 

@m_r2 spring @m_ml @m_m2 m_K m_Z 
@m_r2 spring @m_m2 @m_m0 m_K m_Z 

# Inputs and outputs 
@inl poslnput 0. 

Soutl posOutput @m_m2 

# Add plucking interaction 

Spick nlPluck Sinl §m_ml nlK nlScale 
Listing 3: MIMS description for the model presented in Fig. [ 7 ] 

The graphical UI version of MIMS also provides basic tools 
for generating certain categories of physical structures (i.e., strings, 
membranes, etc.) and performing modal analysis of linear structures. 



Figure 4: The MIMS model editor prototype. 


import("stdfaust.lib"); 
import("mi.lib"); 

inPos = hslider("pos",1,-1,1,0.0001) : si. 

smoo; 

OutGain = 10.; 

m_K = 0.1; 
m_Z = 0.001; 
nlK = 0.05; 
nlScale = 0.01; 


model = ( 

RoutingLinkToMass: 


ground(0.), 


mass (1 .,0 . , 

0 

mass (1.,0 . , 

0 

mass (1 .,0 . , 

0 

poslnput (0 . 

) 


RoutingMassToLink : 


spring (0.05,0.01, 0. 

, 0.) , 

spring(m_K,m_Z, 0., 

0.) , 

spring(m_K,m_Z, 0., 

0.) , 

spring(m_K,m_Z, 0., 

0.) , 

nlPluck(nlK,nlScale) 

t 

par (i, 1,_) 

par(i, 10, _): par(i, 10,! 

), par(i 


with { 

RoutingLinkToMass(10_fl,10_f2,ll_fl, ll_f2, 
12_f1,12_f2,13_f1,13_f2,14_f1,14_f2,ini) 
= 10_fl, 10_f2+ll_fl+13_f2, ll_f2+12_f1+ 
14_f2, l2_f2+13_f1, 14_f1, ini; 
RoutingMassToLink(mO,ml,m2,m3,m4) = mO, ml, 
ml, m2, m2, m3, m3, ml, m4, m2,m3; 

} ; 


process = inPos : model: *(OutGain); 

Listing 4: MIMS description for the model presented in F7gn?e[T] 


4. EXAMPLES AND EVALUATION 

The basic mijaust package contains several examples of virtual in¬ 
struments and use-cases of mass-interaction physics in FAUST. All 
of these examples can be compiled and executed directly as web ap¬ 
plications via the FAUST online editor^jwith generic user interfaces. 
They can also be found as pre-compiled web-apps on the mijaust 
project web-page[^] 

• IPlayTheTriangle: the demonstration model discussed 
previously in this paper (FigureJT]). 

• PolyTriangle: the same model (with a direct force im¬ 
pulse applied instead of a pluck system), using Faust’s abil¬ 
ity to automatically handle polyphonic voice allocation for 
MIDI instruments. 


The FAUST code generated from the model in Code Listing[3]is 
presented in Code Listing [4] The only hand-written element is the 
inPos function, that adds a graphical slider to control the position 
of the input mass. The control-rate output of the slider is smoothed 
to avoid artifacts. 


• PluckedHarmonics: a 150-mass string terminated by two 
fixed points. The first position input allows plucking the 

“https://faust.grame.fr/tools/editor/ 
f https://faust.grame.fr/community/ 
made-with-faust/mi-faust 
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string, and three others are used to press down lightly on the 
string at specific areas in order to bring out natural harmonics. 

• BowedString: a bowed string, using the nIBow interaction. 
The user can control bow pressure and velocity, as well as the 
stiffness of the string. 

• LargeTriangleMesh: a big triangular mesh, fixed at one 
summit, excited by a plucking system and damped by user 
input. 

• Resonator: the audio input is fed into one end of a resonat¬ 
ing physical model. The user can alter the properties of the 
resonator. 

• PhysicalLFO: Using a physical model with slow dynamics 
as a control variable for another synthesis process. Here, the 
wave propagation observed along a very loose string is used 
to modulate the amplitude of a white noise source, generating 
AM modulation going from complex patterns at the onset to 
quasi-sinusoidal modulation as the higher modes of the string 
decay. 

In addition to these examples, two large structures (a 20 by 
30 mass mesh: 20x30mesh and a 1000 mass string: 1000 
massString) were created for model complexity tests. The bench 
test results in Table [T] show the compile time and CPU load for var¬ 
ious models. Large routing functions result in slower compilation, 
and maximum complexity is reached for approx. 1800 physical el¬ 
ements. Overall, fairly complex models run well, with a reasonable 
CPU load. 


5. FUTURE WORKS 

5.1. Faust 

FAUST proves to be well adapted to implement mass-interaction 
physical models. The combination of connection matrices and of 
the use of the letrec environment expression allowed us to seam¬ 
lessly implement the various elements of mi . lib. However, this 
raised some issues that will need to be solved in the future. They are 
presented below. 

5.1.1. Specifying Initial States in letrec 

The letrec environment expression doesn't allow us to specify an 
initial state (i.e., the value of y(n — 1), y(n — 2), etc. at n = 0). We 
got around this problem by implementing the initState function 
which requires some unneeded computation. Hence, letrec could 
be modified to allow this type of expression to be written (rewriting 
Code Listing[X]l : 

equation = x 
letrec{ 

x' = xO; 
x' ' = xl; 

'x = A*x +B*x' + * (C); 

}; 

We believe that this would significantly reduce computation for 
large scale models. 


5. 1.2. Optimizing Routing Matrices 

The current “bare bone” implementation of connection matrices 
(e.g., RoutingLinkToMass in Code Listing [3]» is hard to solve 
by the Faust compiler, preventing large models to be generated (see 
This could be solved by turning this operation into a primitive 
of the language. Compilation time would be significantly reduced 
since pattern matching 0 wouldn’t be involved to solve this type of 
expression. 

6. CONCLUSIONS 

In this paper, we have presented early results of formal integration 
of ID mass-interaction physical modeling into the FAUST environ¬ 
ment, resulting in a new library. The MIMS and physics2faust tools 
allow to automatically generate Faust dsp code for complex topo¬ 
logical models, by expliciting the routing scheme for the model’s 
position and force signals. Several basic models have been imple¬ 
mented and benchmarked, showing promising results. Furthermore, 
Faust’s capabilities offer an efficient solution for playing several 
dynamically allocated and parameter-mapped instances of a physical 
model across large ranges. More generally, this work extends beyond 
mass-interaction modeling and explores the possibilities for describ¬ 
ing complex feedback networks and state space-models in FAUST. 
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Model Name 

N. Masses 

N. Springs 

Faust Comp. Dur. 

CPU Load 

lOOOmassString 

1000 

1002 

- 

- 

20x30mesh 
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1151 

20.576s 

45% 

BowedString 

150 

152 

1.962s 

14% 

IP1ayTheTriangle 

3 

5 

0.029s 

1% 

LargeTriangleMesh 

324 

901 

12.083s 

48% 

PhysicalLFO 

10 

12 

0.032s 

1% 

PluckedHarmonics 

150 

152 

2.192s 

14% 

PolyTriangle 

3 

5 

0.027s 

1% 

Resonator 

30 

32 

0.056s 

4% 


Table 1: Number of masses and springs, compilation duration, and CPU load of the examples. Measurements were made on a Lenovo ThinkPad 
XI Carbon with the following configuration: Linux Manjaro, Intel i7-7500U 4 cores at2.7GHz, 16GiB of RAM, sampling rate of 48KHz, buffet- 
size of256 samples. Programs were compiled as ALSA applications with a GTK interface using faust2alsa. 
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